Generation of spin current by Coulomb drag 
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Coulomb drag between two quantum wires is exponentially sensitive to the mismatch of their 
electronic densities. The application of a magnetic field can compensate this mismatch for electrons 
of opposite spin directions in different wires. The resulting enhanced momentum transfer leads to 
the conversion of the charge current in the active wire to the spin current in the passive wire. 

PACS numbers; 71.10.Pm, 73.63.Nm 



A set of unusual transport phenomena in which 
electron-electron interactions induce transfer of momen- 
tum between distinguishable systems of fermions is 
known as Coulomb drag effect. Conventional Coulomb 
drag [l[ occurs between two spatially separated conduc- 
tors. In the standard setup, see Fig.[Tl dc current Ii flows 
through the active conductor 1 inducing a voltage drop 
V2 in the passive conductor 2. Quantitatively, the effect 
is characterized by the dimensionless drag resistance 



Rd^ lim {e'/h)V2/Ii. 



(1) 



Unlike the usual two-terminal resistance, is sensitive 
to electronic correlations within the conductors. There- 
fore, Coulomb drag effect provides an important tool to 
probe these correlations. Coulomb drag was observed 
experimentally in two-dimensional bilayers [21 and, more 
recently, in one-dimensional quantum wires [3|. 




FIG. 1: Equivalent circuit for measurement of Coulomb drag 
between two quantum wires. Coulomb drag manifests itself 
in the appearence of the potential difference V2 between the 
ends of the open circuit of which the passive wire 2 is a part 
(V2 is positive if it has the polarity indicated). 

A different Coulomb drag- type effect, the spin drag, 
originates in momentum transfer between spin-up and 
spin-down electrons within the same conductor The 
spin drag provides a non-dissipative mechanism of relax- 
ation of a pure spin current. Interactions are therefore 
destructive for spin currents. Because robust generation 
of spin currents is important in view of possible applica- 
tions in spintronics [5[ , the limitations arising due to the 
spin drag effect are now a subject of active research 

In this paper we demonstrate that interactions can in- 
duce spin current rather than suppress it. This is possi- 
ble in a novel type of Coulomb drag effect, interaction- 
induced transfer of momentum between spin-up and spin- 
down electrons that belong to separate conductors. We 



show that this effect can be realized in the standard set- 
ting of Coulomb drag between two clean quantum wires 
in a magnetic field 31 . While the electric current I2 in 
the passive wire is zero, the spin current hs ^ — hi 
can flow [3], i.e. the system acts as a charge current to 
spin current converter. The efficiency of the conversion 
can be characterized by the ratio 



C = l2slh- 



(2) 



Below we show that the drag resistance Rd has a maxi- 
mum at a certain value of Zeeman energy. For 



max{r, |B-Bo|}<Bo 



(3) 



the conversion efficiency C ^ Rd [see Eqs. p8|) and 
([23]) ]. and the dependence of Rd on temperature T is 
described by a power law with the exponent depending 
on the interaction strength, see Eq. (fT5|) . For sufficiently 
strong interaction the power-law dependence crosses over 
to i?d ~ 1 at very low temperatures. We start with a 
heuristic explanation of the origin of the effect, and then 
proceed with the derivation of the results. 

If the electronic densities in the wires ni and 712 were 
equal, the dominant contribution to Rd at low tempera- 
tures would come from processes with large momentum 
transfer between the wires (backscatteringV which may 
result in a finite Rd in the limit T ^ [a H • In re- 
ality, however, the densities are always slightly different, 

\ni - ^2! <C n, n = (rii -|- 712 )/2 

(let us assume that ni < 712), so that the corresponding 
Fermi momenta fci_2 = "'711^2/2 are different as well. In 
this case, the backscattering contribution to Rd is expo- 
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nentially suppressed at low temperatures 

The suppression is easy to understand as follows. To 
the lowest order in the strength of the interwire interac- 
tion, the backscattering contribution to Rd can be writ- 

13, m 



ten as 



JdqJ^^^dioe'-^^HSf^iq^co) 



(4) 



Here L is the length of the region in which the wires 
interact with each other (see Fig. [T]), U2k is 2fc-Fourier 



2 



component of the interwire interaction potential (with 
= (fci + fc2)/2 = W2), and 5f (g,u;) = S,{q,Lo)\^^^^ 
is the Fourier transform of the dynamic structure factor 
Si{x,t) = (pi{x,t)pi{0,0)') (here pi is the local density 
operator for wire i). 
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FIG. 2: (a) Regions in (aj,g)-plane where 51,2 > at T = 
and q ~ 2k. The dark triangle indicates the region where 
5152 > 0. (b) In a magnetic field, the low-energy sectors in 
Si{q,iL!) split in two, which leads to the decrease of To, the 
minimal energy at which Si and S2 overlap at T = 0. 

At T = and q ~ 2k, the two structure factors overlap 
only at lo > Tq v\ki — fc2|, where v = 7rn/2r7i is the 
"average" Fermi velocity, see Fig. [2ja). Because of the 
factor e~^l'^ in Eq. ([4]), this translates to the activational 
temperature dependence of the drag resistance, Rd oc 
g-Tb/T Although at any T > the structure factors are 
finite for all uj and g, the "leakage" of the spectral weight 
beyond the boundaries indicated in Fig. [2l^a) affects only 
the power-law prefactor in the expression for Rd- 

With the backscattering contribution exponentially 
suppressed, Rd is dominated by small momentum trans- 
fer and vanishes at T ^ as i?d cx [l3]. In principle, 
the densities can be fine-tuned to be equal, which would 
increase the backscattering contribution. Another possi- 
bility, which leads to spin current generation, is to place 
the system in a magnetic field. 

In a field the single-particle energies ffeo- of the spin- 
up (t) and spin-down (|) electrons (labeled by cr = ±1) 
include Zeeman contribution b^ka = aB /2. As a result, 
riii > fii^, and the Fermi momenta are 



ki(j — ki — (j6k/2 



(5) 



with Sk{B) ^ B/v (see below). For each wire, the low- 
energy sector in Sf'^{q,Lu) then splits in two, located at 
q — 2fcicr, see Fig. [UJb). The scale Tq is S-dependent 
and vanishes at a certain field Bq, Tq{B) ~ \B — Bq\ (see 
Eq. in]) below). At \B-Bq\ < T, the backscattering con- 
tribution to Rd is no longer exponentially suppressed and 
dominates at sufficiently low temperatures. Moreover, in 
the regime (jSj the main contribution to the integral in 
Eq. ([4]) comes from the overlap of Sii and 82^ In other 
words, almost all of the momentum is transferred from 
spin-down electrons in the active wire to spin-up elec- 
trons in the passive one. Therefore, both Rd and C will 
have a maximum at B — Bq. 

We evaluate Rd and C in the regime ([3]) using the 
bosonization technique [ij]. At energies well below Bq, 
which in turn is small compared with the Fermi energy 



ep, the wire i (i = 1, 2) is described by the Hamiltonian 

+ 9m{^a:'^^Jnf] , (6) 

where m = c, s labels the charge (spin) modes, and the 
bosonic fields satisfy 

[^Ptmix), i9t'7n'{y)] = {i/'2)Sii'S„un' sgn(x - y) . (7) 

For simplicity, we assume that both wires are described 
by the same set of parameters {vm,gm}- These parame- 
ters are related to each other according to 

gc^v/v„ gs{Bo) = l + [2ln{eF/Bo)]-^ (8) 

(so that 1 ~ Qc ^ ffs — 1 > for i?o ^ ef), and the 
velocities Vc> v and Vg < v can be further expressed in 
terms of the interaction within the wires IJ] . 

Fermion operators in the bosonic representation are 

^^Ux) = A^.c«.VPi^e"['''°-(^)+^-^l. (9) 

Here a ~ +1(— 1) for the right (left) moving fermions, 
MiQcr = lAaa ^rc real (Majorana) fermions that sat- 
isfy { 

Mmcn Mi'a'cr' } — ^^ii' ^aa' ^tjtj' (thcsc Operators en 
force correct anticommutation relations between different 
fermionic species) , po ^ Bq jv is the high- momentum cut- 
off, and T]iaa is a linear combination of ipimi'&inu which 
in the leading order in Sq/ef ^ 1 is given by [l5| 

Viaa = V7r/2 [ipic + a'dic + cFLpis + aadig) . (10) 

Fermi momenta /cjg. in Eq. ^ are given by Eq. ([5]) with 
5k{B) = g^B/v,, and Tq{B) (see Fig.^]) at B ^ Bo is 

Tq{B) ^ g,\B - Bq\, BQ^v,\k2-ki\ (11) 

[Bq is the root of the equation gs{B)B — u^lfe — fci|)- 

With the help of Eq. ([9]) , the 2fc-harmonic of the den- 
sity operator pf^ — J^aPfa is written as 



2k 



PoPia- exp [i V27r (i^ic + o-ip^s) + 2ikiax] + H.c. 



where — Pi,-i.api,+i,a- Since the Hamiltonian ([6]) is 
quadratiCj^evaluation of the structure factor is straight- 
forward [14] and yields Si{x, t) = Sia{x, t) with 



S.ia{x,t) = 2pocos(2fcio.a::)]^ 



T/{2pQVm) 

sinhCTrTr^m) 



1 9m/2 



where Tam = x/vm — a{t — iO). 

As discussed above, the condition ^ ensures that the 
main contribution to the integral in Eq. (jj]) comes from 
the nonvanishing overlap of Sii and 5*21; the remaining 
contributions are suppressed as oc exp(— Bq/T). In order 
to evaluate Rd, it is convenient to convert Eq. (|4]) to 
space-time representation, 

= (7r/2)C/|fe / dxdt{it)Si{x,t)S2{x,t). (12) 
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Substituting here 5*11 for 5*1 and and S2] for S'2, we find 



Bo 



\B-Bn 



Bn 



where — C^2fc/27rw and g — {gc + 5s)/2. The function 
F{z) in Eq. (fT^ is given by 



F{z) 



(z/2)3-49 exp(2izC/7r) d^dC, 
n™[cosh(^e + C)cosh(i;-C-C)]' 



26""^, K Z < Zo , (14) 

Z » Zo 

where zg = gc(7''/2)(ws/'Uc) tan[(7r/2)(u5/t;c)] (so that 
zq ^ (1 — 3c)~^ S> 1 for weak interaction). In deriving 
Eqs. (1121) , HH) we changed the integration variables in 
(ini) to ^ = ttTx/vs and C = ttT^, shifted the path of in- 
tegration over C off the real axis by — i7r/2, and evaluated 
the resulting integral in the saddle-point approximation. 
According to Eqs. (fT3|l . (fT4|) . and in agreement with the 
discussion above, Rd{B) has a narrow peak of the width 
5B r-. T <^ Bo aX B = Bf). Its height is given by 



max{i?rf(B)} ^ nXlkL{BoltF){T/Bo) 



43-3 



(15) 



Note that the difference between Vs and Vc is impor- 
tant only at large \B — Bo\ ^T/{1 — gc). In the oppo- 
site limit one can set Vg/v^ — > 1, which yields F{z) — 
\T{g + iz/2TT)f /T'^{2g), in agreement with Eq. the 
corresponding T-dependence is exactly the same as that 
for the drag between two spinless wires [ll|, [lB| • 

In order to relate the conversion efficiency ^ to the 
drag resistance ([T5|) . we note that as far as the passive 
wire is concerned, in the regime ([3]) Coulomb drag in- 
duces the electric field that couples to spin-up electrons 
only. The effect of this field can be described by adding 
to the Hamiltonian of the passive wire a term 



within the leads and that the contacts between the leads 
and the wires are reflectionless; the corresponding con- 
ductance is not affected by the interactions [l3|- The 
condition of vanishing of the total electric current, I2 = 
Iv + he = 0, then yields = -2V^2/5c- Eqs. {I]),© 
and (|17p now give 



C - hs/h - 2{gs/gc)Rd- 



(18) 



Thus, under the conditions ([3]) the dependence of conver- 
sion efficiency C on _B and T is indeed the same as that 
of the drag resistance Rd, as asserted above. 

Eq. (fTS]) does not account for the reduction of Is due 
to the momentum transfer between the two spin subsys- 
tems within the passive wire (spin drag). Indeed, in the 
framework of the Tomonaga-Luttinger model ([6]) the only 
source of spin drag is the backscattering in the spin sec- 
tor, which at T <C -B is exponentially suppressed. The 
dominant contribution to spin drag then comes from the 
processes with small momentum transfer. Accounting 
for these processes requires explicit consideration of the 
nonlinearity of the electronic spectrum [l^ . Proceeding 
along the lines of [l3| , we found the corresponding correc- 
tion to the spin current hs at T <C -B and in the lowest 
non-vanishing order in the interaction strength, 



Shs/hs - -nLo{l - gcfiB/epfiT/Bf 



(19) 



In writing Eq. (jl9|) we took into account that Fermi 
velocities for spin-up and spin-down electrons differ by 
Sv ^ B/k <^ V. The correction (fT9| is small and does 
not affect the validity of Eq. (fT8|) . 

The above consideration is based on the perturbative 
expression Eq. (|4]). In order to analyze the relevance of 
the higher-order contributions, we introduce new fields 



'1/2 



V2, 



:), cf)s = 2~^/'^{ipis+'P2s), 



and similarly defined 9c and Os ■ The fields obey the com- 
mutation relations analogous to Eq. fT]), and their dy- 
namics is governed by the Hamiltonian H = JdxTi. with 



6H2 = e j dx'^d{x)p2]{x) = e J dx ^ {p2c + P2s) , (16) H = [gm^l^x^m)^ +gmidx9m)'^] " 2vXQ{dx(, 



where ^d{x) is drag-induced potential, and p2c and p2s 
are charge and spin densities. The potential ^d{x) 
changes within the region of the length L in which the 
wires interact with each other. Assuming that the wires 
are long. La ^ L, the charge and spin currents in re- 
sponse to 6H2 can be written as 17 1 



AT:vX2kPo cos|\/47r ( 



(l)s) + 2Koxy (20) 



he = (2e7/i) gc (5$<j/2, hs - (2e7/i) 5, 6<^d/2, (17) 

where S^d — ^di—00) — ^d{oo). In writing Eq. pT|) we 
took into account the renormalization of the correspond- 
ing conductances by interactions within the wire [17| . 

On the other hand, the electrostatic potential differ- 
ence V2 induces charge current ly = {2e'^/h)V2- Here 
we assumed that the interactions are efficiently screened 



The second and the third terms here describe, respec- 
tively, the forward and backward scattering between the 
spin-up electrons in wire 2 and the spin-down electrons 
in wire 1, with Xq defined similarly to X2k in Eq. (jl3p . 
andXo-To(B)/i;,. 

The forward scattering term in Eq. (I^D)) leads to small 
corrections to Vc and g^ Sgc/gc ~ —5vc/vc ~ 2g^Ao <C 1, 
which modify the exponent in Eqs. (|13p -(|15 p . 5 — *■ 5 + 
Sgc/2. The backscattering, however, can be relevant in 
the renormalization group sense [l9l |. For L —> 00 and 
Kq ^ it then results in the opening of a gap 



(21) 
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in the excitation spectrum. The gapped state is the 
"zigzag" -ordered state formed by the spin-down electrons 
in wire 1 and the spin- up electrons in wire 2. 

The gap remains open for finite Kq as long as the en- 
ergy gained due to its formation is sufficient to overcome 
the cost of the adjustment of the densities needed to form 
the zigzag order. In the context of quantum wires such 
adjustment (known as commensurate-incommensurate 
transition) was discussed recently in 11,23|- The adjust- 
ment takes place at not too large Kq, Kq < ~ A/v, 
and occurs even when L is finite. As a result, the width 
SB of the peak in Rd[B) saturates at low temperatures, 



5B - max{T, A}. 



(22) 



For L ^ Tj/A the zigzag order can not be formed and 
Eq. (fT5|) is applicable. In this case max{i?t;(i?)} <C 1 for 
all T . The higher-order contributions become important 
for L > w/A and at r < A d, i, [13, While finding 
the detailed dependence Rd{T) in this regime is beyond 
the scope of this Letter, the limiting values of Rd and C 
at T ^ can be found as follows. Imagine that the two 
wires are connected to noninteracting reservoirs and a 
bias is applied only to the electrons with spin a in wire i. 
The resulting current of electrons with spin a' in wire i' 
is /i'a' = Giia'^irjVia, whcrc Gi'a'^icr = Gia,i'a' IS the Corre- 
sponding conductance. At T — > the spin-up electrons in 
wire 2 are "locked" with the spin-down electrons in wire 
1, and we expect that Gix.ii, G2|,2Tj ^11,21 ^ e^/2/i. At 
the same time, Gii^ij, G'2x,2i —> e^//i, while Gi^,2\ 0. 
Setting Via = Vi, = I, 



HI, 



> e /h, while G'i|.2X 
we find 



Rd ^ 1/4, C ^ 1/2. 



(23) 



To conclude, we showed that in the presence of the ap- 
plied magnetic field the standard Coulomb drag measure- 
ment setup acts as a charge current to spin current con- 
verter. Both the drag resistance and the conversion effi- 
ciency exhibit a maximum at a certain value of the field 
controlled by the density mismatch between the wires. 

Our results are applicable for long {kLo 1) ballis- 
tic quantum wires. The wires studied in fs^ exhibit a 
well-defined conductance quantization, which guarantees 
that the elastic mean free path exceeds the length of the 
wires Lq. While it is very plausible that kLo ^ 1 for at 
least some of the samples studied in ^3] (with Lq ranging 
from 0.4 /im to 4/im), the density of electrons in these 
wires is difficult to estimate. Fortunately, such estimate 
is available for the coupled-wire system studied in 
L w Lq ~ 10 /xm and kLo ~ 10'^. Although the exper- 
iments [2l'| focus on the momentum-resolved tunneling, 
the same system can be employed to study the Coulomb 
drag effect as well. For this system, the typical density 
mismatch |ni — n2\/n ~ 10~^ corresponds to Bo ~ IK 
(which amounts to the applied field of ^ 3 Tesla) , hence 
the regime ([3]) is well within the reach of the experiments. 
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